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1 Introduction 

On the one hand, partial differential-integral inclusions (PDIIs in short) play 
an important role in characterizing many social, physical, biological and en- 
gineering problems, one can see Balasubramaniam and Ntouyas [3] and refer- 
ences therein. On the other hand, reflected backward stochastic differential 
equations (RBSDEs in short) associated to a multivalued maximal monotone 
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operator defined by the subdifferential of a convex function has first been in- 
troduced by Gegout-Petit [B]. Further, Pardoux and Racanu [18J gave the 
existence and uniqueness of the solution of RBSDEs, on a random (possi- 
bly infinite) time interval, involving a subdifferential operator was proved 
in order to to give the probabilistic interpretation for the viscosity solution 
of some parabolic and elliptic variational inequalities. Following, Ouknine 
[H], N'Zi and Ouknine [13j, Bahlali et al. [U, [2] discussed this type RBS- 
DEs driven by Brownian motion or the combination of Brownian motion and 
Poisson random measure under the conditions of Lipschitz, locally Lipschitz 
or some monotone conditions on the coefficients. El Karoui et al. [S] have 
got another type RBSDEs different here, where one of the components of the 
solution is forced to stay above a given barrier, which provided a probabilis- 
tic formula for the viscosity solution of an obstacle problem for a parabolic 
PDF. Since then, there were many works on this topic. One can see Matoussi 
[12], Hamedene [HIIS], Lepeltier and Xu [H], Ren et al. [121120] and so on. 

The main tool in the theory of BSDEs is the martingale representation 
theorem, which is well known for martingale which adapted to the filtration 
of the Brownian motion or that of Poisson point process (Pardoux and Peng 
|17j . Tang and Li [25]). Recently, Nualart and Schoutens [H] gave a martin- 
gale representation theorem associated to Levy process. Furthermore, they 
showed the existence and uniqueness of solutions to BSDEs driven by Teugels 
martingales associated with Levy process with moments of all orders in ^15] . 
The results were important from a pure mathematical point of view as well 
as in the world of finance. It could be used for the purpose of option pricing 
in a Levy market and related partial differential equation which provided an 
analogue of the famous Black-Scholes partial differential equation. 

Motivated by the above works, the purpose of the present paper is to 
consider RBSDEs associated to the subdifferential operator of a lower semi- 
continuous convex function driven by Teugels martingales associated with 
Levy process which considered in Nualart and Schoutens [HIITS]. We obtain 
the existence and uniqueness of the solutions for RBSDEs. As its application, 
we give a probabilistic interpretation for a class of PDIIs. 

The paper is organized as follows. In Section 2, we introduce some prelim- 
inaries and notations. Section 3 is devoted to the proof of the existence and 
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uniqueness of the solutions to RBSDEs driven by Levy processes by means of 
the penahzation methods. A probabihstic interpretation for a class of PDIIs 
by our RBSDEs is given in the last section. 

2 Preliminaries and notations 

Let T > be a fixed terminal time and JF, P) be a complete probability 
space. Let {Lt : t G [0,T]} be a M- valued Levy process corresponding to 
a standard Levy measure u whose characteristic function has the following 
form: 

E{e^^^^) = exp[iaut — ^a'^uH + t J^{e'''^^ — 1 — iuxl[\x\<:i})iy{dx)], 

where a G M, a > 0. Furthermore, Levy measure u satisfying the following 
conditions: 

(1) U^^yXdy) <oo- 

(2) j]_££[c e'^^y^i'(dy) < oo, for every 5 > and for some A > 0. 

This shows that Lt has moments of all orders. We denote by (i/^*))j>i 
the linear combination of so called Teugels Martingale y/*"* defined as follows 
associated with the Levy process {Lt : t G [0,T]}. More precisely 

J^t — -\-Ci^i^iIt -\- ■ ■ ■ -\- Ci^ilt ) 

where y}'^ = Lf - E\l}p] = Lf - tE[Lf] for alH > 1 and Lf are power- 
jump processes. That is, = Lt and Lf' = ^^^^^^(ALt)* for i > 2. It 
was shown in Naulart and Schoutens [S] that the coefficients Ci^k correspond 
to the orthonormalization of the polynomials 1, x^, ■ ■ ■ with respect to the 
measure fj.{dx) = x'^v{dx) + a'^dQ^dx): 

Qi—l C-i iX -\- Ci i—\X ~l~ ■ ■ ■ ~l~ Cj^x. 

The martingale {H^^^)i>i can be chosen to be pairwise strongly orthonormal 
martingale. Furthermore, [H^''\ H^^'^],i ^ j, and {[H^^\ H^^^]t — t}t>o are 
uniformly integrable martingales with initial value 0, i.e. < H^^\H(j) >t= 
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Remark 2.1. If v = 0. we are in the classic Brownian case and all non-zero 
degree polynomials qi{x) will vanish, giving = 0,^ = 2,3,-- - . ///i only 
has mass at 1, we are in the Poisson case; here also H^^^ = 0, i = 2, 3, • • • . 
Both cases are degenerate in this Levy framework. 

Let jV denote the totality of P— null sets of J^. For each t e [0,r], we 
define 



Let us introduce some spaces: 

•H? = {cp : .Fj— progressively measurable process, real-valued process, 
s.t. E /q \ipt\^dt < cxo} and denote by the subspace of formed by the 
predictable processes; 

•iS^ = : jFj— progressively measurable process, real-valued process, s.t. 
^(supo<t<rlv5Wn < oo}; 

•P — {{xn)n>o :be real valued sequences such that Yl°loxf < oo}. We 
shall denote by H'^{P) and V'^{f') the corresponding spaces of /^-valued pro- 
cess equipped with the norm 



Now, we give the following assumptions: 
(HI) The terminal value ^ e L'^{Q,Tt, P); 

(H2) The coefficient / : [0,T] x Q x Rxf R is .Ft -progressively 
measurable, such that /(•,0,0) e H^; 

(H3) There exists a constant C > such that for every {cu,t) e Q x 
[0,T],(|/i,zi), (2/2,^2) eRxP 

\f{t,y,,zr) f{t,y„z,)\' < C{\y, - y,\' + \\z, - z^W'); 

(H4) Let $ : i? ^ (— oo,-|-oo] be a proper lower semicontinuous convex 
function. 

Define: 

Dom{^) := {x e R: < +00}, 

d^{x) := {x* e R :< X* , X - u > +^{u) < ^{x),\/x G R}, 
Dom{d^)) ■.= {xeR:d^^ 0}, 

Gr{d^)) := {ix,x*) e R^ : x e Dom{d^)),x* e d^x)}, 



Tt = cr{L„0< s<t)Vj\f. 
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Now, we introduce a multivalued maximal monotone operator on R de- 
fined by the sub differential of the above function $. The details appeared in 
Brezis jl]. 

For all X G i?, define 

= miuj^d |x - + $(?/)). 

Let J„(x) be the unique solution of the differential inclusion x G Jn{x) + 
^9$(J„(x)), which is called the resolvent of the monotone operator A = 5$. 
Then, we have 

Proposition 2.1. (1) '■ R ^ R is Lipschitz continuous; 

(2) The Yosida approximation of is defined by An{x) := V$„(a;) = 
n{x — Jn{x)),x G R which is monotonic and Lipschitz continuous and there 
exists a G interior{Dom{^)) and positive numbers R,C satisfies 

< V$„(2;)*, z-a>> R\An{z)\ - C\z\ - C, for all z e R and n E N; (1) 

(3) For all x E R, An{x) E A{Jn{x)). 

Further, we assure that ^ E Dom{^) and E^{^) < oo. 
This paper is mainly discuss the following reflected backward stochastic 
differential equations. In doing so, we first give its definition. 

Definition 2.1. By definition a solution associated with the above assump- 
tions (^, /, $) is a triple (Yt, Zt, Kt)Q<t<T of progressively measurable pro- 
cesses such that 

Yt = i + /(s, Z,)ds + KT-Kt- ET=i j7 zf^dHf\ 0<t<T, 

(2) 

and satisfying 

(1) {Yt,Zt)o<t<T E X V^iP) and {Yt,0 < t < T} is cddldg (right 
continuous with left limits) and take values in Dom{^); 

(2) {Kt,0 < t < T} is absolutely continuous, Kq = 0, and for all pro- 
gressively measurable and right continuous process {{at,(3t),0 < t < T} E 
Gr{d^), we have 

^Yt-at){dKt + fitdt) < 0. 
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In order to obtain the existence and uniqueness of the solutions to (2), 
we consider the following penalization form of (2): 

yr = i + j^[f{s.Y:,z:)-A^{Y:)]ds~YZJlzf^^''dH'i\Q<t<T, 

(3) 

where ^, / satisfies the assumptions stated above and An is the Yosida ap- 
proximation of the operator A = 5$. Since An is Lipschitz continuous, it is 
known from the result of [TBJ, that Eq. (3) has a unique solution. 

Set K]^ = — < t < T. Our aim is to prove the sequence 

(y", Z^, i^") convergence to a sequence {Y, Z, K) which is our desired solu- 
tion to RBSDEs (2). 

3 Existence and uniqueness of the solutions 

The principal result of the paper is the following theorem. 

Theorem 3.1. Assume that assumptions on and $ hold, the RBSDEs 
(2) has a unique solution (Yf, Zf, Kt)o<t<T- Moreover, 

lim^^oo ^ supo<t<T \yr - ^t? = 0, 

lim™E \\Z^ - ZtW^dt = 0, 
lim^^oo supo<t<T \Kt - Kt\^ = 0, 
where {Y^,Z^) he the solution of Eq. (3). 

In the sequel, C > denotes a constant which can change from line to 
line. 

The proof of the Theorem is divided into the following Lemmas. 

Lemma 3.1. Under the assumptions of Theorem 3.1, there exists a constant 
Ci > such that for all n> 1 

i?(supo<,<r \Yr\' + \Z^\'ds + lAniYDlds) < Ci. 

Proof. By Ito formula, we have 

\Yr-a\' = |e-a|2 + 2/f(n"-a)/(s,y,",Z,")c?s 

-2/f (y« - a)An{Y,-)ds - Ji'WZ^W'ds (4) 

-2E.=i/f(n"-a)z«'"rf//:. 
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Taking expectation on the both sides and considering (1), we obtain 
E\Y^^-a\'^ + E j^WZl'W'^ds 

< m-a\' + 2Ej^{Y--a)f{s,Y-,Z-)ds 
-2REJ^ \Ar,(YJ')\ds + 2CEJ^ \Yp\ds + 2C. 



(5) 



Further, we get 

- ap + ^ J^WZ^W^ds + 2RE \AniY^'')\ds 

< E\^-a\' + 2Ej^^{Y--a)f{s,Yr,Z-)ds 
+2CEj^\Yp\ds + 2C (6) 

< E\^ - ap + CE \Y^ - aprfs + \E WZ^W'^ds + C, 

where we have used the elementary inequahty 2ah < j3'^a^ + ^ for all a, 6 > 0. 
So, we have 

E\Y^ - a\'^ + ]^E WZ^W^ds <C{1 + E lY^"" - a\^ds). (7) 
Gronwall inequality shows 

^^l^t" - < C, Vn. (8) 

So that 

E|yj"|2 < C, Vn. (9) 
Prom (6) and (7), it is easy to show 

E r iWZ^W + \A^{Y:^)\)ds < C, Vn. (10) 

^0 

Bulkholder-Davis-Gundy inequality shows the desired result of Lemma 3.1. 

□ 

Lemma 3.2. Under the assumptions of Theorem 3.1, there exists a constant 
C2 > such that for all n> 1 

Ej^\A^{Yr)\'ds<C2. 
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Proof. Without loss of generality, we assume $ > 0, $(0) = 0. Let tl^n — 
For the convexity of and Ito formula, we have 



Taking expectation on the both sides, we obtain 



= E^^iO + E VMYDfis, Y:\ Z^)ds 
-^E \A^{Yr)\'ds. 



(11) 



(12) 



Using the elementary inequality 2ab < Pa^ + ^, we obtain 
EMY,-) + '-E \A^{Yr)\'ds 

< EMO + \A^{Y-)\^ds + ^E J^^ \f{s,Y;,Z:)\'ds. (13) 

Further, we have 

EUYn + 'tiE lAr^iYDl'ds 

< EMO + ^EJ^ lA^YDl^ds + ^EJ^ \f{s,Yr, Z-)\^ds 

< EMO + hE^I \AniY-)\''ds + |/(.,n", Z':)\Hs 

< EMO + i-Ej^ |A„(F«)prf. + f E/f m^ds 
+^Ej^'\\Z:\\'ds + ^Ej^\f{s,OM'ds. 



(14) 



Lemma 3.1 shows 

EMYr) + ^„Ej^\Ar,{Yr)\'ds<^, 
which implies the desired result. □ 

Lemma 3.3. Under the assumptions of Theorem 3.1, {Y"',Z'^) be a Cauchy 
sequence in S"^ x V^{f). 

Proof. By Ito formula, we have 

|yn _ ym|2 ^ jT | _ Z^W^ds 

= 2 j^iYr - Yr)[fis, y:\ z^) - f{s, rr, zT)]ds 

-2 ff{Y^ - yr)(^n(n") - A^{Yn)ds (15) 

-2 E^=l /t^(n" - i;'")(^i'^'" - zP'"')dHi. 



For 

we can obtain 

-{Yr - YDiA^Y^ds - AUYD) < i^\An{Yr)\' + l\AUYn\'- 
So, we get 



< 2CE - + l^s" - YrWZ": - ZT\)ds 

< 2Ci?/f [(1 + - Yr\' + ^^\Z- - ZT\']ds 
+EL^{£\An{Yrr + i^\AUYr)\')ds. 



(16) 



Choosing P such that ^ < |, then 

supo<,<ri?|yr - yrP + i^^/o" ll^r - ^riP^^^ < + 

Further, Bulkholder-Davis-Gundy inequahty shows 

^ sup |y"-y™P + iE /" < c(- + — ), (17) 

o<t<T 2 io 

which shows our desired result. □ 

In order to obtain the existence of solutions to Eq. (2), we give the 
following Lemma appeared in Saisho |21j . 

Lemma 3.4. Let {K^"'\ n E N} be a family of continuous functions of finite 
variation on . Assume that: 

(i) sup„ < a < oo, < t < oo; 

(ii) lim„^oo i^^") = K E C([0, +oo); R). 

Then K is of finite variation. Moreover, if {f^'^\n E N} is a family of 
continuous functions such that lim„^oo /'■"■' = / ^ ^([Oj R)i then 
lim„_oo £ < fL^'KdK:^ >= < A, dKu>, yO<s<t< oo. 

Proof of Theorem 3.1 Existence. Lemma 3.3 shows that (!"", Z") is a 
Cauchy sequence in space S"^ x V'^ip). We denote its limit as (F, Z). At the 
same time, it is easy to verify that {K^)neN converges uniformly in L^(i7) to 
the process K. = lim„_^oo /g An{YJ^)ds, that is i?supo<t<T \Kt ~ ^tl"^ = 0- 

Denote H^{0, T; R) be the Sobolev space consisting of all absolutely con- 
tinuous functions with derivative in L^(0,r). Lemma 3.2 shows 
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SUPneAf -^11 -^"111^2 (o,T;R) < 

which imphes that the sequence is bounded in if^(0, T; i?)). So, 

there exists an absolutely contimious function K G LP'{VL; H^{0, T; R)) which 
is the weak limit of K^. Further, ^ = Vt, where -Vt G d^{Yt). 

In the following, we verify that {Y, Z, K) is the unique solution to Eq. 
(2). Taking a subsequence, if necessary, we can show 

SUPo<t<Tl^r-^t| ^0, 

^'^Vo<t<T\K't - Kt \ -^0. 

It follows that Yt is caglag and Kt is continuous. Further, let {ex, {3) be a 
caglag process with values in Gr(9$), it holds 

< J„(y(") - au dK^ + (3tdt >< 0. 

Since J„ is a contraction and uniformly converges to Y on [0, T]. It follows 
that J„(y^") converges to pr{Y) uniformly, where pr denotes the projection 
on Dom($). Lemma 3.4 shows 

< pr{Yt) - at, dKt + >< 0. 
In order to complete the proof of the existence, we need to verify the following 

P{Yt e Dom{^), < t < T) = 1. 
For the right continuity of y, it suffices to prove 

P{Yt e Dom{^)) = 1, < i < r. 

If not so, there exists < t < T and Bq ^ T such that P{Bq) > and 
Yt^ ^ Dom{^)yuj e Bq. By the right continuity, there exists (5 > 0, Bi e 
such that P(Si) > 0, Yt ^ Dom{^) for (cj, t) e Bi x [to, to + 5]. 

Using the fact ^^gn ^ lo \^n0^p)\ds < oo and Fatou Lemma, we get 

/^^ /;;+'liminf„^oo \An{Yp)\dsdP < oo, 

which is impossible for liminf„_>oo \An{YJ^) | = oo on the set Bi x [to, + 

Uniqueness. Let (Yt, Zt, Kt)o<t<T and (F/, Z[, K^)o<t<T be two solutions 
for Eq. (2). Define 

{AYt, AZt, AKt)o<t<T = {Yt - F/, Zt - Z[, Kt - K[)o<t<T. 
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Ito formula shows 

H 



E|Ayi|2 + £;/; WAZtW'ds 



= 2E /; AYs[f{s,Y,,Z,)-f{sX,K)]ds + 2E // AYJAK^. (18) 

Since 9$ is monotone and G d^{Yt), G d^(Y^'), we obtain 

^ Ay.dA/sT, < 0. 

Further, 

E\AYt\'^ + E WAZtW^ds < CEJ^ {AY^l'^ds + ^E WAZ^W^ds. 
Gronwall inequality shows the uniqueness of the solutions to Eq. (2). 

4 Applications 

In this section, we study the link between RBSDEs driven Levy processes 
and the solution of a class of PDIIs. Suppose that our Levy process has 
the form of Lf = at + Xt where Xt is a pure jump Levy process with Levy 
measure u{dx). 

In order to attain our main result, we give a Lemma appeared in [T5] . 



Lemma 4.1. Let c : x [0, T] x M — > M 5e a measurable function such that 

|c(s,?/)| < a,(2/2 A \y\) a.s., 

where {a^, s G [0, T]} is a non-negative predictable process such that E a^ds < 
oo. Then, for each < t < T, we have 

+ !l k^i.s^y)^i.dy)ds. 
Consider the following coupled RBSDEs: 

Yt = h{LT) + f{s, Ls, Zs)ds + KT-Kt- !l z'i^dHf\ 0<t<T, 

(19) 

where E\h{LT)\'^ < oo. 
Define 

u^it, X, y) = u{t, x + y)- u{t, x) - |^(t, x)y, 
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where u is the solution of the following PDIIs: 

'^{t,x) + f{t,u{t,x), 

+ J^u\t,x,y)u{dy) e d^u{t,x)), (20) 



f (i, x) + a'^it, x) + fit, u(t, x), {u^\t, 



u{T,x) = h{x) e Dom{^), 
where a' = a + /{|j^|>i} yi^{dy) and 

^^'\t, x) = u\t, X, y)pi(yHdy) + x)(J^ y'i^(dy))'^, 
and for i >2 

u^^(t,x) = J^u\t,x,y)pi(y)u(dy). 

Suppose that u is C^'^ function such that ^ and |^ is bounded by poly- 
nomial function of x, uniformly in t. Then we have the following 

Theorem 4.1. The unique adapted solution of (19) is given by 
Yt^u{t,Lt), 

Kt = J^Vsds, -Vs e d^u{t,Lt)), 

= lR^J'^it^Lt-,y)pi{y)iy{dy), i > 2, 
= J^u\t,L,_,y)pM^{dy) + 

Proof. Applying Ito formula to w(s, L^), we obtain 
m(T, Lt) - u{t, Lt) 



(21) 



Lemma 4.1 apphed to n(s, Lg- + y) — u{s, LgJ) — ||(s, Ls-)y shows 

= E»=i (/m ^*^(s, Ls-,y)pi{y)u{dy))dH'i^ + 1^1(5, Ls-,y)u{dy)ds. 

(22) 

Note that 

= f/') + tL;Li = ( / y^iy{dy))-^Hi'^ + tEL^, (23) 

where ^Li = a + J{^y^>^ V^idy). 

Hence, substituting (22) and (23) into (21) yields 
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h{LT)-u{t,Lt) 

+ f^[u\s, L,_,y)p,{y)v{dy) + f (s, y^v{dy))^dH^''^ 

+ E»=2 (/r«H^> y)Pi{y)y{dy))dH'i\ 
Prom which we get the desired result of the Theorem. □ 
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